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Casimir Forces and Geometry in Nanosystems 

Thorsten Emig^^ 

Casimir interactions, predicted in 1948 [15, 14] between atoms and macroscopic 
surfaces, and probed in a series of high precision experiments over the past decade 
[46, 50, 8], are particular important at micro-meter to nano-meter length scales due 
to their strong power-law increase at short separations between particles. Therefore, 
in constructing and operating devices at these length scales it is important to have an 
accurate understanding of material, shape and geometry dependence of these forces. 
In particular, the observation of Casimir forces in devices on sub-micron scales has 
generated great current interest in exploring the role of these forces for the development 
and optimization of micro- and nano-electromechanical systems [9, 17, 16]. These 
systems can serve as on-chip fully integrated sensors and actuators with a growing 
number of applications. It was pointed out that Casimir forces can make an important 
contribution to the principal cause of malfunctions of these devices in form of stiction 
that results in permanent adhesion of nearby surface elements [10]. This initiated 
interest in repulsive Casimir forces by modifying material properties as well as the 
geometry of the interacting components [11, 43, 52]. 

The study of fluctuation induced forces has a long history. When these forces re- 
sult from fluctuations of charges and currents inside particles or macroscopic objects, 
they are usually summarized under the general term van der Waals forces [55]. This 
interaction appears at the atomic scale in the guise of Keesom, Debye, London, and 
Casimir-Polder forces. An important property of all this interactions is their non- 
additivity: The total interaction of macroscopic objects is in general not given by 
the sum of the interactions between aU pairs of particles that form the objects. This 
inherent many-body character of the force leads to interesting and often unexpected 
behavior but renders the study of these forces also a difficult problem. Commonly 
used approximations as pairwise additivity assumptions become unreliable for systems 
of condensed atoms. The collective interaction of condensed macroscopic systems is 
better formulated in terms of their dielectric properties. Such a formulation was es- 
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tablished by Lifshitz for two parallel and planar, infinitely extended dielectric surfaces 
[47], extending Casimir's original work for perfect metals. In practice, one encoun- 
ters objects of finite size with curved surfaces and/or edges, like structured surfaces, 
spheres or cylinders. Also, in small-scale devices often more than two objects are at 
close separation and one would like to know the collective effects resulting from non- 
additivity. In this chapter we shall encounter a selection of examples for interesting 
behavior of fluctuation forces that result from shape and material properties that have 
been obtained from a recently developed method that makes it possible to compute van 
der Waals-Casimir interactions for arbitrary compact objects based on their scattering 
properties for electromagnetic waves [32, 33, 30, 42]. 

1.1 

Casimir effect 

The Casimir effect is the attraction between two uncharged, parallel and perfectly 
conducting plates [14]. For this simple geometry the interaction can be obtained 
directly from the plate induced change of the energies of the quantum mechanical 
harmonics oscillators associated with the normal modes of the electromagnetic field. 
The derivation given here follows closely the one originally presented by Casimir. 
Consider two parallel and planar surfaces of size L x L and separation d. We assume 
that the system is at zero temperature so that the interaction is given by the ground state 
energies of harmonic oscillators. When we are interested in the pressure (force per 
plate area L^) between large plates with L ^ d,we can ignore edge effects ~ L and 
allow for a continuum of wave vectors parallel to the plates. For a perfect conductor the 
tangential electric field has to vanish at the surface and the normal modes correspond 
to the allowed wave vectors k = (k||,7rn/rf) where k|| is the two-dimensional wave 
vector parallel to the plates . The hnear dispersion of photons yields the eigenfrequencies 
Wnfcii = c^fcy + {■nn/dY so that the ground state energy becomes 

^ = / rf'^^ll 2'^"^ll ' (1-1) 

n=0 ^ ^ 

where we have included a factor of 2 since for each mode with n ^ exist two 
polarizations. The primed summation assigns a weight of 1/2 to the term for n = 
0. Obviously, the expression of Eq. (1.1) is divergent. This is a consequence of 
the assumption that the surfaces behave as a perfect conductor for arbitrarily high 
frequencies. In practice, as pointed out by Casimir, for very high frequencies (X-rays, 
e.g.) the plates are hardly an obstacle for electromagnetic waves and therefore the 
ground state energy of these modes will not be changed by the presence of the plates. 
We implement this observation by introducing a cut-off function x{z) that is regular at 
z = Q with x(0) = 1 and vanishes, along with all its derivatives, for 2 ^ oo sufficiently 
fast. After a change of variables, ui = cJk'i + {-Kn/dy with c^A;||dA;|| = Lodio, we 



1.1 Casimir effect 3 



obtain the finite expression for the energy 

S=— ^/(n) with /(n)=/ uj\{^/uJc)cLj (1.2) 

where i^c is a cut-off frequency. As mentioned before, we are interested in the change of 
the energy due the presence of the plates. Let us imagine that we increase the separation 
d between the plates to infinity, thus creating empty space. When we subtract the energy 
of the latter configuration from the total energy of Eq. (1.2), we obtain the change in 
energy that is the relevant interaction potential between the plates. When the separation 
d tends to infinity, the sum in Eq. (1.2) can be replaced by an integral, yielding after 
the substitution Q, = ■Knc/d the energy 



h 
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L^d / dnf{Q,) with /(Q) = / w^xWwc)dw ■ (1.3) 
Jo Jn 



As expected, the energy E^o is proportional to the volume L'^d of empty space and to a 
cut-off dependent factor that is given by the integrals of Eq. (1.3). This factor describes 
the self-energy of the bounding surfaces. It is infinite for perfect conductors which 
correspond to u)c — > oo. For a non-ideal conductor or any other material this factor is 
finite but depends on material properties hke the plasma wavelength for a metal. Now 
we compute the change in energy when the plates are moved in from infinity, 



E /(-) - / 



dnf{n) 



(1.4) 



The difference between the sum and the integral is given by the Euler-Maclaurin 
formula Er=o7(n) - dn}{n) = -^^/'(O) + ^/'"(O) + O(/"(0)). This series 
of derivatives of odd order can be truncated in the Umit of perfect conductors oic — > oo 
since /'(O) = 0, /'"(O) = -2{'Kc/df and/'^^CO) ~ (c/d)^(c/da;c)""^ TheCasinur 
potential hence becomes 



and the pressure for perfect metal plates is 



_ _Tr^ he 
1?~ 240 # 



(1.6) 



The interesting fact is that the amplitude of the interaction is universal, i.e., independent 
of the cut-off that can be viewed as a simplified description of a real metal. This implies 
that for any pair of surfaces with metallic response in the limit of small frequencies 
a; ^ the interaction at asymptotically large separations is described by the potential 
ofEq.(1.5). At a separation of d = 100nmEq.(l. 6) yieldsapressixre of 1.28-10""' atm 
or 13.00 Pa. 
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1.2 

Dependence on Shape and Geometry 

Casimir interactions result from a modification of the fluctuation spectrum of the 
electromagnetic field due to boundaries or coupling to matter. This suggest that these 
interactions strongly depend on the shape of the interacting objects and geometry, 
i.e., relative position and orientation. The most commonly encountered geometry is a 
sphere-plate setup that was used in the first high-precision tests of the Casimir effect [46, 
50]. This geometry has been successfully used ever since in most of the experimental 
studies of Casimir forces between metalhc surfaces [60, 27, 17, 21, 19, 25, 20, 51]. In 
order to keep the deviations from two parallel plates sufficiently small, a sphere with a 
radius much larger than the surface distance has been used. The effect of curvature has 
been accounted for by the "proximity force approximation" (PFA) [55]. This scheme 
is assumed to describe the interaction for sufficiently small ratios of radius of curvature 
to distance. However, this an uncontrolled assumption since PFA becomes exact only 
for infinitesimal separations, and corrections to PFA are generally unknown. 

At the other extreme, the interaction between a planar surface and an object that is 
either very small or at an asymptotically large distance is governed by the Casimir- 
Polder potential that was derived for the case of an atom and a perfectly conducting plane 
[15]. There have been attempts to go beyond the two extreme limits of asymptotically 
large and small separations by measuring the Casimir force between a sphere and a 
plane over a larger range of ratios of sphere radius to distance [44]. 

Until very recently, no practical tools were available to compute the electromagnetic 
Casimir interaction between objects of arbitrary shape at all distances, including the 
important sphere-plate geometry. Progress in understanding the geometry dependence 
of fluctuation forces was hampered by the lack of methods that are applicable over a 
wide range of separations. Unlike the case of parallel plates, the eigenvalues of the 
Helmholtz equation in more complicated geometries are in general not known and a 
summation over normal modes, as in the original Casimir calculation of Section 1 . 1 , is 
not practical. Conceptually, the effects of geometry and shape are difficult to study due 
to the non-additivity of fluctuation forces. 

For some decades, there has been considerable interest in the theory of Casimir forces 
between objects with curved surfaces. Two types of approaches have been pursued. 
Attempts to compute the force explicitly in particular geometries and efforts to develop 
a general framework which yields the interaction in terms of characteristics of the 
objects like polarizability or cixrvature. Within the second type of approach, Balian 
and Duplantier studied the electromagnetic Casimir interaction between compact and 
perfect metallic shapes in terms of a multiple reflection expansion and derived also 
explicit results to leading order at asymptoticaUy large separations [4, 5]. For parallel 
and partially transmitting plates a connection to scattering theory has been established 
which yields the Casimir interaction of the plates as a determinant of a diagonal matrix 
of reflection amplitudes [40]. For non-planar, deformed plates, a general representation 
of the Casimir energy as a functional determinant of a matrix that describes reflections 
at the surfaces and free propagation between them has been developed in Ref. [31]. 
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Later on, an equivalent representation has been applied to perturbative computations in 

the case of rough and corrugated plates with finite conductivity [45, 54, 53]. 

Functional determinant formulas have been used also for open geometries that do 
not fall into the class of parallel plates with deformations. For the Casimir interaction 
between planar plates and cylinders, a partial wave expansion of the functional deter- 
minant has been employed [36, 6]. More recently, a new method based on a multipole 
expansion of fluctuating currents inside the objects has been developed [32, 33, 30]. 
This method allows for accurate and efficient calculations of Casimir forces and torques 
between compact objects of arbitrary shape and material composition in terms of the 
scattering matrices of the individual objects. A similar scattering approach has been 
developed in Ref . [42] . 

In this Section three examples for the strong geometry dependence of Casimir forces 
will be made explicit. First, an overview on forces between deformed and structured 
surfaces will be given. The interactions are obtained from both a perturbative and 
numerical evaluation of a functional determinant representation of Casimir interactions 
between ideal metal surfaces. As a second example we describe the effects that occur 
in the interaction of one-dimensional structures as cylinders and wires and related 
non-additivity phenomena for more than two objects. Finally, as an example for the 
interaction between compact objects, the Casimir force between metallic spheres is 
presented for the full range of separations, covering the crossover from the asymptotic 
Casimir-Polder law to proximity approximations. The analysis of the last two examples 
is based on a scattering approach. 

1.2.1 

Deformed surfaces 

The dependence of the Casimir force on shape and material properties offers the 
opportunity to manipulate this interaction in a controlled way, e.g., by imprinting 
patterns on the interacting surfaces. It has been shown that a promising route to this 
end is via modifications of the parallel plate geometry [38, 34, 35]. The corrections due 
to deformations, such as sinusoidal corrugations, of the metal plates can be significant. 
In searching for non-trivial shape dependences, Roy and Mohideen [61] measured 
the force between a sphere with large radius and a sinusoidally corrugated plate with 
amplitude a « 60nm and wavelength A « l.l/xm. Over the range of separations H « 
0.1 — 0.9 /im, the observed force showed clear deviations from the dependence expected 
on the basis of decomposing the Casimir force to a sum of pairwise contributions (in 
effect, an average over the variations in separations). Motivated by this experiment, the 
effect of corrugations on the Casimir force between surfaces has been studied without 
using pairwise additivity approximations. The analysis is based on a path integral 
quantization of the fluctuating field with appropriate boundary conditions which leads 
to a functional determinant representation of the Casimir energy [38J which can be 
evaluated perturbatively for a small deformation amplitude [34, 35] or numerically for 
general amphtudes [28, 12, 13]. In recent experiments [18] the Casimir force between 
a gold sphere and a silicon sixrface with an array of nanoscale, rectangular corrugations 
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has been measured and the results were found to be consistent with the theory based 
on a numerical evaluation of functional determinants for ideal metals (see below). 
Qualitative agreement can be expected only when material properties are taken into 
account in addition to shape. 




Fig. 1.1 Configuration of a flat and a corrugated plate at mean separation H. 

While we will be interested in the interaction between flat and corrugated surfaces 
as depicted in Figs. 1.1, 1.4, first of all we consider two surfaces with arbitrary uni- 
axial deformations without overhangs so that their profiles can be described by height 
functions ha{yi) (q = 1,2 for the two surfaces), with J dyi ha{yi) = 0. It is fur- 
ther assumed that the surfaces are perfectly conducting and infinitely extended along 
the plane spanned by yy = (5/1,3/2). As explained before for two planar plates, the 
Casimir energy at zero temperature corresponds to the difference of the ground state 
energies of the quantized electromagnetic field for plates at distance H and at 00, 
respectively. To obtain this energy, we employ a path integral quantization method. For 
general, non-uniaxial deformations or objects of more general shape, it is necessary to 
consider the action for the electromagnetic field since the two polarizations (TM for 
transverse magnetic waves and TE for transverse electric waves) are coupled. How- 
ever, for the uniaxial deformations under consideration here, we can develop a simpler 
quantization scheme, by a similar reasoning as used in the context of waveguides with 
constant cross-sectional shape [34]. In this case, the two polarizations are independent 
modes which do not couple under scattering between the surfaces. For TM waves all 
field components are then fully specified by a scalar field corresponding to the electric 
field along the invariant direction, 

$TM(t, 3/1,2/2,2:) = E2{t,yi,y2,z), (1.7) 

with the Dirichlet boundary condition c&tm | Sq = on each surface Sa ■ The TE waves 
are analogously described by the scalar field 

*TE(t,2/i,y2,2) = B2{t,yi,y2,z), (1.8) 

with the Neumann boundary condition dn^TE\sa = 0, where dn is the normal deriva- 
tive of the surface Sa pointing into the space between the two plates. After a Wick 
rotation to the imaginary time variable X" — id, both fields <1?tm and <1>te can be 
quantized using the Euclidean action 

S{^} = IJ d^X{\/<i>f . (1.9) 
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In order to obtain the change in the ground state energy that is associated with the 

presence of the plates, we now consider the partition functions Zo and 2n for the scalar 
field Euclidean action both with Dirichlet (D) and Neumann (N) boundary conditions 
at the surfaces. We implement the boundary conditions on the surfaces Sa using delta 
functions, which leads to the partition functions 

1 r ^ 

Zu = — V^llll5mXc)]eM-S{^}/'i}, (1-10) 

1 r ^ 

2n = — V<S>llllS[d^<^>{X^)]eM-Sm/K) , (1.11) 

where Zq is the partition function of the space without plates. Here Xi (y) = [y , hi {yi )] 
andX2(y) = [y, H + (2/1)], where y = (2/0,2/1,2/2) = (2/0, y|i), and 2/0 = ict, is a 
parametrization of the plates in 4D Euclidean space. The Casimir energy £ per unit area 
(at zero temperature) that results from moving the plates in from infinity is obtained 
from the partition function as 

£{H) = E{H) - lim E{H) , (1.12) 

with 

fir 

E{H) = -—l\nZn + lnZj,], (1.13) 

where A is the surface area of the plates and the limit where the overall Euclidean length 
in time direction, L, tends to infinity is implicitly assumed. The partition functions can 
be expressed as fimctional determinants, using auxiliary fields (for details see Ref. [35]), 

In^^D = - ^ IndetMD, ln.2N = - ^ IndetMn . (1.14) 
The kernels Md and Mn are given by 

[MD]„/3(y,y') = [ffa(2/i)]'/'G[X„(y)-XHy')][sMyi)]'^', (1-15) 
[MN]a/3(y,y') = [ffa(2/i)]'/'5„„(„)9„^(,.)G[X„(y)-X^(y')] 

xMyi)Y^\ (1.16) 

where ga{yi) = 1 + [h'a{yi)]'^ is the determinant of the induced metric, and ndyi) = 
(— l)"5a ^^^(2/i)[/ia(2/i), 0, —1] is the normal vector to the surface Sa, while 

G(x) = ^^ (1.17) 

is the free Euclidean space Green's function with x = (y, z). Eqs. (1.12) to (1.17) 
constitute the functional determinant representation of the Casimir interaction. This 
representation is exact. To proceed, the determinant has to be evaluated either by 
perturbation theory in the deformation amplitude or numerically for a specific shape of 
the surface profiles. First, we present the perturbative approach. 
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(1.19) 



For both boundary conditions {X = D, N) we devide by the partition function 2^x,oo 
for ^ oo and expand ln(Zx /^x,oo ) in a series \n{Zx /Zx,oo) \ a+ln{Zx/Zx,oo) \ i + 
ln{Zyi/Zx,oo)\2 + ■ ■ ., where the subscript indicates the corresponding order in ha. 
The lowest order result is 

ln(2x/2x.^)|„ = ^^ (1.18) 

for both types of modes, corresponding to two flat plates, cf. Eq. (1.5). The first order 
result ln(.Ex/-Zx,oo)|i vanishes since we assume, without loss of generality, that the 
mean deformations are zero, / dy\ha{y\) = 0. The second order contribution is given 
by 

ln(2x/2x,oo)|2 = ^ ^ f d'y {[h,{y,)f + [^(t/i)]'} 
-Ijd'yj dPy'Kx{y-y') \j^[hi{y,) - h,{y\)f + \[h2{y,) - h2{y\)f 
-IJd^J A'Qx(y-y') [hi{yi)h2(y[) + h2{yi)hi{y[)] . 

The terms in the first row are local contributions which are identical for TM and TE 
modes. They follow also from a pairwise summation approximation (PWS) that sums 
a "renormalized" Casimir-Polder potential over the volumes of the interacting bodies 
[35]. The remaining terms are non-local and cannot be obtained in approximative 
schemes. For Dirichlet boundary conditions, the kernels depend only on |y — y'| and 
are given by 

. s 1 TT^ 1 COsh^(s) 

^ / N 7r^ 1 sinh^(s) ^, 

«-(^) = 128 ^'-''^ 

where s = ■Ky/{2H). The kernels for Neumann boundary conditions assume a more 
complicated form since the normal derivative breaks the equivalence of space and 
time directions. Hence, they depend separately on |yo — 2/ol and lyy — y|||. Their 
explicit form can be found in [35]. The results obtained so far apply to general uniaxial 
deformations of both surfaces. 

Now we apply these results to the important case of corrugated plates. We begin 
with the geometry depicted in Fig. 1.1 which is parametrized by 

hi{y{) = acos(27rj/i/A) , and /i2(2/i) = 0. (1-22) 

For this profile, the computation of the partition function to second order in a reduces 
to the Fourier transforming of the kernels with respect to yi. The corresponding 
expression for £ in Eq. (1.12) can be written as 

£ = £o + £ci, (1.23) 
where £o is the energy per unit area of two flat plates [see Eq. (1.5)] and 



_ hca 



+ C>(a'), (1.24) 
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where the index of of Eci stands for corrugated-flat geometry. The functions that 
describe the A-dependence in this expression can be computed exactly [35]. They can 
be expressed in terms of the polylogarithm function Li„(z) — J^J^^ z" jv^'^ , leading to 



480 
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30 ^ ' 1920a; ^ ^ 24 ^ ' 



XX _L 
;^Li3(u) + — -Li4(K) + Li5(u) + 
24 62tt 04%'' 



2567r3a 



Li6(u) 



TT 

945, 
(1.25) 



TEia;) = 
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48 



+ 



1440 
l + 2x2)Li2C 
7 _ . . . 1 
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Li3(i, 



1287r2 



Lis (it) + 



2567r3a 



(1.26) 



with u = exp(— 47ra;). Figure 1.2 displays separately the contributions from Gtm and 
Gte to the corrugation induced correction £ci to the Casimir energy. While Gtm (H/X) 
is a monotonically increasing function of H/X, Gte{H/X) displays a minimum for 
H/X « 0.3. 




H/?v 

Fig. 1.2 Rescaled correction Scf to the Casimir energy due to the 
corrugation as given by Eq. (1 .24) (upper curve). The lower curves 
show the separate contributions from TM and TE modes. The 
rescaiing of £rj is chosen such that the corresponding prediction 
of the pairwise summation (PWS) approximation [corresponding to 
the iocai terms of Eq. (1.19)] is a constant (dashed iines). 



Examining the limiting behaviors of Eq. (1.24) is instructive. In the limit A » //, 
the functions Gtm and Gte approach constant values, and the total Casimir energy 
takes the A-independent form 

H^ 720 V H^ 



(1.27) 



Note that only in this case both wave types provide the same contribution to the total 
energy and the result agrees with the pairwise summation approximation (see Fig. 1.2). 
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In the opposite limit of A <C -ff both Gtm and Gte grow Unearly in H/X. Therefore, 
in this limit the correction to the Casimir energy decays according to a slower power 
law in H, as 

with an amplitude proportional to 1/A. Note that this behavior is completely missed 
by the pairwise simimation approach which always yields a A independent Casimir 
energy in the presence of modulations on one plate [35]. As we will discuss below in 
the context of the numerical approach, the apparent divergence for A — > in Eq. (1.28) 
is an artifact of the perturbative expansion which assumes that the amplitude a is the 
smallest length scale. 

Next we turn to a numerical approach for computing the functional determinants 
of Eq. (1.14). Such an approach has been developed for periodic surface profiles 
in [28, 12], In this approach, it is convenient to compute directly the Casimir force 
F = —OhS per unit area which is the sum of TM and TE contributions, F = Ftm + ^te 
that according to Eq. (1.14) are given by (for X = D,N) 

Fy, = Tr (Mx^a^Mx) . (1.29) 

The right hand side of this expression is always finite, and no divergences due to self- 
energies have to be subtracted. The trace in Eq. (1.29) can be computed efficiently by 
Fourier transforming M with respect to y, y'. The so transformed operator can then 
be transformed to block-diagonal form by making use of the periodicity of the surface 
profile along the j/i direction. In this representation the blocks can be numbered by 
the wave vector qi £ [0, 27r/A) along the yi direction. A block matrix with label qi 
couples only waves whose momenta differ from the Bloch momentum qi by integer 
multiples of 27r/A. The integers multiplying 27r/A number the matrix elements within 
a block matrix. Hence, the problem of computing the total trace has been simplified 
to the computation of the trace of each block matrix with label qi . Finally, integration 
over qi from to 27r/A and over the unrestricted momenta qo, q2 (along the time 
direction and invariant spatial direction of the surfaces, respectively) yields the force 
of Eq. (1.29). For the particular choice of a rectangular corrugation [see Fig. 1.5(a)] 
analytic expressions for all matrix elements of M can be obtained. For details of the 
implementation of the numerical approach and expressions for the matrix elements see 
[12]. 

In analogy to the profile of Fig. 1 . 1 we consider the corresponding situation of a flat 
plate and a plate with a rectangular corrugation profile parametrized by 

h,{y,) = I , , (1.30) 

' \ -a for A/4< < A/2 

and continuation by periodicity /ii(2/i) = -I- nA) for any integer n. The numerical 
results for the total Casimir force between the two plates is shown in Fig. 1 .3 for different 
corrugation wavelengths A. For all A, the forces at a fixed separation H are bounded 
between a minimal force Fao and a maximal force Fq. For small A/a the upper bound 
Fo is approached whereas for asymptotically large A/o the force converges towards the 
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lower bound Fao ■ Analytic expressions can be derived for these bounds. For large A, the 
corrugated surface is composed of large flat segments with a low density of edges. At 
sufficiently small surface separations <C A the main contribution to the force comes 
irom wavelengths which are much smaller than the scale A of the surface structure. Thus 
in the dominant range of modes diffraction can be neglected, and a simple proximity 
force approximation [55] should be applicable. Such an approximation assumes that 
the total force can be calculated as the sum of local forces between opposite flat and 
parallel small surface elements at their local distance H — hi{yi). No distinction is 
made between TM and TE modes. This procedure is rather simple for the rectangular 
corrugation considered here since the surface has no curvature (except for edges). There 
are only two different distances H + a, H — a which contribute one half each across the 
entire surface area, leading for A ^ oo to the proximity approximation for the force, 



/A — TT ?lC 1 

= / ~ ~ W 2 



1 1 

+ 



(H-a)4 (H + a) 



(1.31) 



In the limit A ^ the important fluctuations should not get into the narrow valleys of 
the corrugated plate. Even for small but finite A this picture should be a good, though 
approximate, description since it still effects the wavelengths of order H which give the 
main contribution to the force. Thus one can expect that the plates feel a force which is 
equal to the force between two flat plates at the reduced distance H — a. Fortunately, 
this expectation can be checked by an explicit calculation since the leading part of 
determinant of Mx in the hmit A ^ can be computed. Indeed, this computation 
confirms the expectation, leading to the Casimir force per surface area [12] 

with equal contributions from TM and TE modes. Notice that this result is not analytic 
ina/H and is exact in the limit A ^ 0. As we have seen before, perturbation theory 
for smoothly deformed surfaces always yields corrections to the interaction of order 
a^. However, for small a/H, the result of Eq. (1.32) has the expansion 



^'1 



(1.33) 



which indicates that perturbation theory is not applicable if A ^ a. This implies that 
the apparent divergent behavior for A ^ in Eq. (1.28) actually disappears for A ~ a. 



1.2.2 

Lateral forces 



As a natural generalization of the geometry of the previous section, we study the Casimir 

interaction between two sinusoidally corrugated plates. For direct correspondence to 
experiments for this type of configuration [22], we consider the specific profiles 

hi{yi) = acos(27ri/i/A) , and h2{yi) = acos{2n{yi + b)/\) , (1.34) 
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Fig. 1.3 Total Casimir force as function of the mean plate 
separation H. Shown is the relative change of the force compared 
to the total Casimir force Ffl^t between two flat plates. The two 
boid curves enclosing the numerical data are the anaiytical resuits 
Fq for a ^ (upper curve) and Fac tor \ ^ oo (lower curve), 
see text. 



Fig. 1.4 Geometry used for calculating the lateral Casimir force 
between two corrugated plates with lateral shift b. The equilibrium 
position is at 6 = A/2. 

which are shifted relative to each other by the length b (see Fig. 1.4). When these 
profiles are substituted into the general expression for the second order term of the 
partition function of Eq. (1.19) one finds for the Casimir energy 



with £ci given in Eq. (1.24), and where the corrugation-corrugation interaction energy 
fee can be calculated in terms of the kernels Qx{y) in Eq. (1.19). Besides oscillating 
contributions to the normal Casimir force from £cc{b), a lateral force 




£ — £() + 2i£ef + £i 



(1.35) 




(1.36) 
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is induced by the corrugation-corrugation interaction. This lateral force is much better 

suited for experimental tests of the influence of deformations, since there is no need 
for subtracting a larger baseline force (the contribution of flat plates) as in the case of 
the normal force. The lateral force can be also employed as a actuation mechanism in 
mechanical oscillators as we will see in Section 1 .4. In analogy to the previous section, 
the corrugation-corrugation interaction can be expressed as 

Jtm ( "T ) + -^TE 



Sec — 



hca 



/ 271-6 



with 



120 



16x* 



l) arctanh(v'it) + \At 



(1.37) 



+ f2*('^'3'^)+lfc^(«'^'^)+3i 



*(w,6,-) 
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24 



+ 6^*^'^'''^)+ 256^ 



$(«,3, 



*(«,6,-) 



32n 



3 X i 

^ +2 + 8to 



Jtm (a;) + Jte{x) 



47r^ 



1. 
(1.38) 



(1.39) 



where w = exp(— 47ra;) and $(0, s, a) = •*'°/('^ + '^)'' '^^ Lerch transcendent. 

In the limit of large corrugation length, HjX 0, this result agrees to lowest order 
with a pairwise summation approximation where Jtm(O) + Jte(O) = 7r^/120. 

At the other extreme of \ <^ H, Jtm (a) -j- JTE(a;) decays exponentially fast. This 
decay distinguishes the lateral force from the normal force. In particular, for large 
X = H/\, we get to leading order 



(1.40) 



Since Jtm (a:) + JTE(a;) is positive for all values of x, the equilibrium position of two 
modulated surfaces is predicted at 6 = A/2. This corresponds to aligning the maxima 
and minima of the two corrugations (cf Fig. 1.4). 

The numerical approach for computing the functional determinant in the case of 
periodic surfaces can be also applied to the lateral force [13]. We consider again a 
rectangular corrugation but now on both surfaces with a lateral shift of b, cf. Fig. 1.5. 
The numerical results for the lateral force in this geometry are summarized in Figs. 1.5 
and 1.6. Fig. 1.5 shows the numerical result for the lateral force for a shift b = A/4 and 
different values of A/a over more than 4 orders of magnitude for the gap 6 = H — 2a, 
together with two approximate results (PFA and PWS) and the perturbative result for 
sinusoidal profiles for A ^ i/. An exponential decay of the force as predicted by 
perturbation theory can be clearly observed. 

The PFA yields a lateral force per unit area Fiat.ppA = [2£o{H) — £o{H — 2a) — 
£o{H + 2a)]/A for < 6 < A/2 where £0 has the same meaning as before. Fiat.PFA 
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(a) ,y^ 



(b) 
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Fig. 1.5 (a) Geometry consisting of two parallel plates with 
laterally shifted uniaxial rectangular corrugations, (b) Lateral force 
Flat (in units of normal force Fq between flat surfaces) atb = A/4 
for the geometry shown in (a) as function of the gap 5 (solid 
curves). Plotted are also the proximity force (PFA, dash-dotted 
curves) and pairwise summation (PWS, dashed curves) 
approximations, and the perturbative result Fpt that follows from a 
calculation for sinusoidal profiles (dotted curves). 

changes sign at b — A/ 2 discontinuously which is an artifact of this approximation. The 
pairwise summation (PWS) of Casimir-Polder potentials is strictly justified for rarefied 
media only but it is often also applied to metals, using the two-body potential U (r) — 
— {■K/24)hc/r' with the amplitude chosen such as to reproduce the correct result for flat 
ideal metal plates [7]. It yields a lateral force fiat, pws = — ^ /y^ rf^x d^x'i/dx — 
x' I ) with Vi and V,- denoting the semi-infinite regions to the left and right of the two 
surfaces in Fig. 1 .5(a), respectively. Fiat.pws can be obtained by numerical integration. 
For small gaps 5, both approximations agree and match the exact numerical results. 
Beyond 5 > A/20 the PFA staits to fail since it does not capture the exponential 
decay of Fiat for increasing 5. The PWS approach has a slightly larger validity range 
and reproduces the exponential decay. However it deviates by at least one order of 
magnitude from Fiat for 5 > 2.5A. 

Although the perturbative result of Eq. (1.37) applies to sinusoidal surfaces, it is 
instructive to compare it to the numerical results for the rectangular profiles. Since the 
lateral force decays exponentially. Fat ^ e~^'"^^^, with the characteristic scale set by 
the modulation wavelength of the profile, the force at large H should be determined 
by the lowest harmonic of the periodic surface profile. This implies a universal lateral 
force for large H ^ \ that is independent of the precise form of the surface coiTugation. 
This universal force is the force between two sinusoidal surfaces where the amplitude 
follows from the projection of an arbitrary periodic profile of wavelength A onto a 
sinusoidal profile with the same wavelength. The latter force follows from Eq. (1.37) 
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Fig. 1 .6 Shape dependence of Fiat on the lateral surface shift b 
at fixed distance H = 10a for different corrugation lengths. The 
dashed and the dotted curves represent the PWS and the full 
perturbative result for sinusoidal profiles with arbitrary H/X, 
respectively. 



and in the limit A <C -ff is given by 



where we assumed an amplitude ao for the sinusoidal profiles. For the rectangular 
corrugation of Fig. 1.5 the lowest harmonic has the amplitude ao = Aa/n. When we 
compare Fpt and the numerical results of Fig. 1 .5(b) we find indeed excellent agreement 
for distances 5 > \. 

The universal behavior of the lateral force is also clearly demonstrated by the de- 
pendence of the lateral force on the surface shift 6. Corresponding numerical results 
together with PWS approximations and the force that follows from the full perturbative 
result of Eq. (1.37) for sinusoidal surfaces with arbitrary H/\ are shown in Fig. 1.6 for 
fixed H = 10a and varying A/a. With decreasing A, three regimes can be identified. 
For A 2> //, the force profile resembles almost the rectangular shape of the surfaces, 
and the PWS approximation yields consistent results. For smaller A, yet larger than 
H, the force profile becomes asymmetric with respect to 6 = A/4 and more peaked, 
signaling the crossover to the universal regime for \ < H where the force profile 
becomes sinusoidal. In the latter case, for not too small A/a ^ 10, the numerical 
results for Fiat indeed agree well with the perturbative result for sinusoidal surfaces 
with arbitrary H/X. We note that the PWS approach fails to predict the asymmetry of 
the force profile, and the PFA even predicts no variation with 6 for < 6 < A/2. To 
observe this universal behavior of the lateral force experimentally, one should consider 
surfaces with very small corrugation wavelengths in the range of nanometers so that 
the exponential decay does not diminish the force for // ^ A too strongly. 




(1.41) 
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1.2.3 

Cylinders 

In this Section we give examples for two central aspects of fluctuation forces: Ef- 
fects resulting from the non-additivity and the particular properties of systems with a 
codimension of two which plays a special role as we will see below. These problems 
are considered in the context of interactions between cylinders and sidewalls. It has 
been demonstrated that Casimir forces in these geometries have only a weak logarith- 
mic dependence on the cylinder radius [36] and can be non-monotonic [59, 57, 56] 
- consequences of codimension and non-addivity. These forces between quasi-one- 
dimensional structures could be probed in mechanical oscillators that are composed of 
nano-wires or carbon nanotubes. Exact results for the interaction can be obtained by 
employing a recently developed scattering approach for Casimir forces [32, 30] . This 
approach is based on the concept that electromagnetic Casimir interactions result from 
fluctuating currents inside the bodies. It is possible to formulate an effective action for 
the multipole moments Qa ,x of the currents inside the bodies where a labels the bodies 
and X is a multi-index that numbers polarizations (electric and magnetic multipoles) 
and the elements of the basis for the multipole expansion, e.g., cyhndrical waves. The 
effective action can then be written as the quadratic form 

S = X! X! Q;,xM„<,,,xx'Qa',x' , (1.42) 

a, a' X,X' 

with the matrix kernel 

M«„-,xx' =«:{[(Ta)"']^^,5„<,/ -U„„.,xx'(l-<5««')} , (1-43) 

where k is the Wick rotated frequency, ui = icK, the matrix is the so-called T- 
matrix of object a that relates incoming and scattered waves and Vaa' is a "translation" 
matrix that relates the incoming wave at object a to the outgoing wave at object a' . 
The T-matrix is related to the scattering matrix of the object, §„, by the relation 
Ta = {Sa — l)/2. Analytic results for all elements of the scattering matrix are 
available for symmetric shapes such as cylinders and spheres. The matrix contains 
all information about shape and material composition of the object that is relevant to the 
Casimir interaction. The translation matrices Uq,q,/ are independent of the properties of 
the interacting bodies and depend only on the relative position (separation vector) of the 
objects a and a' and the properties of the fluctuating field. For the electromagnetic field, 
the translation matrices are known in many bases, e.g., for cylindrical and spherical 
waves [30], To obtain the Casimir energy, the multipole fluctuations are integrated out, 
leading to the determinant of the infinite dimensional matrix M. Integration over all 
frequencies k yields the interaction energy 

^ he [°° , , detM 

£ = 7^ dnln - — , (1.44) 

where the division by the determinant of the matrix Moo accounts for the subtraction of 
the residual energy of the configuration where the separations between all objects tend 
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to infinity. Since the translation matrices decay to zero with increasing separation, the 

matrix Moc is given by Eq. (1 .43) with the Vaa' set to zero. In the special case of two 
objects, the energy can be simplified to [32] 



he 
2tt 



/ dKlndet(l-N), (1.45) 
Jo 



where N = T1U12T2U21. 

First, the scattering approach is applied to two parallel, infinitely long, perfectly 
conducting cylinders of equal radius R and center-to-center separation d, see Fig. 1.7. 
For this geometry it is most convenient to use cyUndrical vector waves for the multipole 
expansion. This basis consists of the vector fields M^^°^(x) = x V''-°-'(x) 
for magnetic (M) multipoles and N'^'°^(x) = x V x V''°'(x) for electric 

(E) multipoles where q — ^ {(^/c)^ — and incoming (i) and outgoing (o) waves 
differ in the definition of the vector fields V'(x) = zJ,„(gr)e™"*e*'=-^ V°(x) = 
zHm\qr)e'"^'^ e''''''^ . Here {r,(f),z) denote cylindrical coordinates and Jm, are 
Bessel and Hankel functions of the first kind. In this basis the matrices of Eq. (1.43) 
assume a simple form where the multi-index X represent now the polarization (M or E), 
the wave vector fcz along the cylinder axis and the partial wave index m. The T-matrix 
is diagonal in polarization, kz and m with diagonal elements 

where we have applied a Wick rotation lj — icK which leads to modified Bessel 
functions of the first (/„) and second (Kn) kind. The translation matrices are diagonal 
in polarization and kz, and the elements are identical for both polarizations, 

2 _ 

Ul2,Mkznm = U\2,Ekznm = — (" *) Km.-n{pd) (1-48) 

ITT 

2 _ 

f^21, Mfc-nm = t^21, Bfc.nm = — * Km-nipd) (1-49) 

with p = \/k^ -\- k'j. Due to the decoupling of electric and magnetic multipoles, 

corresponding to transverse magnetic (TM) and transverse electric (TE) field modes, 
respectively, the Casimir energy of Eq. (1.45) has two independent contributions, 
£2 = £tm + £te, with 

'Stm{te) = ^ y pdplndet(l - Ntm(tb)) , (1-50) 

where L {-^ 00) is the cylinder length and the determinant runs here only over the 
partial wave indices m, m! of the matrix elements 

T.^f^--ipd)^^K^,,.i,d) (1.51) 
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This result for the energy can be also obtained from a scalar field theory where TM 
(TE) modes correspond to Dirichlet (Neumann) boundary conditions [56, 36]. 

At large separations R, only matrix elements with m = m' = for TM modes 
and m — m! = 0, ±1 for TE modes contribute to the energy. When the determinant 
in Eq. (1.50) is restricted to these elements, we find for the interaction of two cylinders 
for large d/ Rio leading order 



^™ ~^'^hdnJ{d/R)) V ~ InW^ + 



(1.53) 



The asymptotic interaction is dominated by the contribution from TM modes that 
vanishes for ii ^ only logarithmically. 



2R 




(d-2R)/R 



Fig. 1.7 Casimir energy for two cylinders of equai radius /? as a 
function of surface-to-surface distance d- 2R (normalized by the 
radius). The energy is divided by the PFA estimate 
'^pfa''^' = -jll^hcL^R/{d - 2iJ)5 for the energy which is 
applicable in the limit d 2R only. The solid curves show 
numerical results; the dashed lines represent the asymptotic 
results of Eq. (1 .53). The inverse logarithmic correction to the 
leading order result for TM modes cause very slow convergence. 
For the parameter range shown here it was sufficient to consider 
m = 40 partial waves to obtain convergence. 



For arbitrary separations higher order partial waves have to be considered. The 
number of partial waves has to be increased with decreasing separation. A numerical 
evaluation of the determinant and integration has revealed an exponentially fast conver- 
gence of the energy in the truncation order for the partial waves, leading to the results 
shown in Fig. 1.7 [56]. It should be noted that the minimum in the curve for the total 
electromagnetic energy results from the scaling by the PFA estimate of the energy. The 
total energy is monotonic and the force attractive at all separations. 
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The interaction between cylinders is very distinct from the Casimir or van der Waals 
interaction which is reported in literature [55]. Usually, the interaction is proportional 
to the volumes of the interacting objects, i.e., for two spheres of radius R where the 
Casimir energy ~ B!' /(f . This scaling with volumes follows also from a pairwise 
summation of two-body forces. However, from the interaction of two parallel plates 
one knows that the interaction can scale also with the surface area. This two examples 
would suggest for two parallel cylinders of length L an interaction energy ~ LR^/d^' 
or ~ LR^ /d!^. But the actual results of Eq. (1.53) has a much weaker, only logarithmic 
dependence on the radius. Is is interesting to look at the variation of the decay exponent 
of d for the Casimir energy as a function of the codimension of the object. The exponent 
is (—3, — 2 + e, —7) for codimensions 1 (plates), 2 (cylinders), 3 (spheres), respectively, 
and hence not monotonic. For a codimension of two, the Casimir interaction is most 
long ranged. The physical reason for the unexpected scaling of the cylinder interaction 
is explained by considering spontaneous charge fluctuations. On a sphere, the positive 
and negative charges can be separated by at most distances of order R <^ d. The 
retarded van der Waals interactions between the dipoles on the spheres lead to the 
Casimir-Polder interaction [15]. In the cylinder, fluctuations of charge along the axis 
of the cylinder can create arbitrary large positively (or negatively) charged regions. 
The retarded interaction of these charges (not dipoles) gives the dominant term of the 
Casimir force. This interpretation is consistent with the difference between the two 
types of polarizations, since for TE modes such charge modulations cannot occur due 
to the absence of an electric field along the cylinder axis, cf. Eq. (1.53) and Fig. 1.7. 
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Fig. 1 .8 Casimir force between two cylinders parallel to one plate 
or sandwiched between two plates vs. the ratio of sidewall 
separation to cylinder radius {H - R)/R, affixed distance 
d = AR between the cylinders, normalized by the total PFA force 
per unit length between two i solated cylinde rs, 
FpFA = -|(?ic7r^/1920)V-ff/(d- 2R)7. The solid lines refer 
to the case with one plate, while dashed lines depict the results for 
two plates. Also shown are the individual TE (blue) and TM (red) 
contributions to the force. 
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As a second example, the effect of sidewalls on the interaction of two cylinders is 
considered. The geometry consisting of either one or two plates at a separation H from 
the two cylinders is shown in Fig. 1.8. For this type of geometry the mean stress tensor 
has been computed numerically and it has been observed that the force between two 
one-dimensional structures changes nonmonotonically when H is increased [59, 57]. 
This many-body effect can be studied by the scattering approach. Instead of studying 
directly the interaction of the cylinders and plates via their T-matrices, it is more 
convenient to employ the method of images to describe the effect of the sidewalls 
[30, 56]. For perfectly conducting sidewalls, their effect on the electromagnetic field 
can be taken into account by replacing the free space Green's function by a half- 
space or slab Green's function. This results in an expression for the Casimir energy 
similar to Eq. (1.44) that depends only on the T-matrices of the two cylinders and 
translation matrices that connect the original cylinders and their mirror images. The 
expression of the energy can be computed again numerically by truncating the partial 
wave expansion at sufficiently high order. The resulting Casimir force between two 
cylinders with one or two sidewalls as function of the sidewall separation H is shown 
in Fig. 1.8. Two interesting features can be observed. First, the attractive total force 
varies non-mono tonically with H: Decreasing for small H and then increasing towards 
the asymptotic limit between two isolated cylinders for large H, cf. Eq. (1.53). The 
extremum for the one-sidewall case occurs a.1 H — 0.27 R, and for the two-sidewall 
case is at iif — i? « 0.46Ji. Second, the total force between the cyUnders for the two- 
sidewall case in the proximity limit H Ris larger than for H/R — > oo. As might 
be expected, the //"-dependence for one sidewall is weaker than for two sidewalls, and 
the effects of the two sidewalls are not additive: not only is the difference from the 
H ^ oo force not doubled for two sidewalls compared to one, but the two curves 
actually intersect. 

A simple generic argument for the non-monotonic sidewall effect has been given in 
[57]. It arises from a competition between the force from TE and TM polarizations as 
demonstrated by the results in Fig. 1.8. An intuitive perspective for the qualitatively 
different behavior of the TE and TM force as a function of the sidewall distance 
is obtained from the method of images. For the TM polarization (corresponding to 
Dirichlet boundary conditions in a scalar field theory) the Green's function is obtained 
by subtracting the contribution from the image so that the image sources have opposite 
signs. Any configuration of fluctuating TM charges on one cylinder is thus screened by 
images, more so as H is decreased, reducing the force on the fluctuating charges of the 
second cylinder This is similar to the effect of a nearby grounded plate on the force 
between two opposite electrostatic charges. Since the reduction in force is present for 
every charge configuration, it is there also for the average over all configurations. 

By contrast, the TE polarization (corresponding to Neumann boundary conditions in 
a scalar field theory) requires image sources of the same sign. The total force between 
fluctuating sources on the cylinders is now larger and increases as the plate separation 
H is reduced. Note, however, that while for each fluctuating source configuration, the 
effect of images is additive, this is not the case for the average over all configurations. 
More precisely, the effect of an image source on the Green's function is not additive 
because of feedback effects: the image currents change the surface current distribution. 
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which changes the image, and so forth. For example, the net effect of the plate on the 
Casimir TE force is not to double the force as H R. The increase is in fact larger 
than two due to the correlated fluctuations. 

A similar but weaker non-monotonic dependence on H of the force between the 
cylinders is also observed for separations d that are different from the particular choice 
in Fig. 1.8. Also, the force between the cylinders and the sidewalls is not monotonic 
in d but the non-monotonicity is then smaller since the effect of a cylinder on the force 
between two bodies is smaller than the effect of an infinite plate. 



1.2.4 
Spheres 



So far geometries with a direction of translational invariance have been considered. 
In the limit of ideal metal surfaces, this invariance lead to a decoupling of the two 
polarizations of the electromagnetic field. Any geometry of experimental interest will 
obviously lack this symmetry beyond some length scale. Hence, it is important to 
study geometries without this sjmimetry. Compact objects of arbitrary shape obviously 
do not have an invariant direction. Therefore, the two polarisations are coupled and 
the matrices in Eq. (1.43) assume a more complicated form. A natural choice for 
a basis are now vector spherical waves for which the translation matrices Vc,a',xx' 
carry an index X — {E or M, I, m) which represents polarisation E or M and the order 
I > 1, m = —I, . . . ,1 of the spherical waves. In contrast to the cylindrical matrices 
of Eq. (1.48), the translation matrix couples now E and M polarisation and all matrix 
elements are known explicitly [30]. 

Here we focus on the simplest case of two compact objects: two perfect metal spheres 
of equal radius R and center-to-center separation d, see Fig. 1.9. The T-matrix of a 
dielectric sphere is known from the Mie theory for scattering of electromagnetic waves 
from spherical particles. Due to spherical symmetry, the E and M polarisations for all /, 
m are decoupled so that the T-matrix is diagonal and the coupling of polarisations occurs 
through the translation matrices only. After a Wick rotation to imaginary frequency 
uj = icK the matrix elements assume in the perfect metal Umit the form 



TufMlml'm' — (~l)'7r"r7 — ~r ^Su'Smm' (1-54) 

7,+ i(k7?) + 2k7?/;^i(k7?) 

TEElml'm' - (-1) ^ J^i^^^^^kR) + 2kRK[1,^ {kR) " ^^'^^^ 

Substitution of these matrix elements together with those of 1Uq,q,' from [30] in Eq. ( 1 .45) 
yields the Casimir energy of two spheres. For asymptotically large d, the energy has 
only contributions from I = I' = 1 (dipoles) and one obtains the Casimir-Polder 
interaction between two polarizable particles [15] where the electric and magnetic 
dipole polarizabilities of a perfect metal sphere are given by = J?^ and um = 
—R^/2. This result can be extended to smaller separations by including higher order 
multipoles with I > 1 that generate higher powers of R/d. One obtains the asymptotic 



/ Casimir Forces and Geometry in Nanosystems 




Fig. 1 .9 Casimir energy of two metal splieres, divided by the PFA 
estimate Spfa = -(7r^/1440)?ic/J/{d - 2Rf, which holds only 
in the limit Fl./d 1/2. The label I denotes the multipole order of 
truncation. The curves i = oo are obtained by extrapolation. The 
Casimir-Polder curve is the leading term of Eq. (1 .56). Inset: 
Convergence with the truncation order I for partial waves at short 
separations. 



series [32] 

„ hcR<^^ fRy ^, _ 

^^--Vd^Z:="(^dJ . (1-56) 

n=0 ^ ^ 

where the first 8 coefficients are co = 143/16, ci=0, C2=7947/160, C3=2065/32, 
C4=27705347/100800, cs = -55251/64, ce = 1373212550401/144506880, ct = 
—7583389/320. The energy at all separations can be obtained by truncating the 
matrix N defined below Eq. (1.45) at a finite multipole order I, and by computing the 
determinant and the integral numerically. The result is shown in Fig. 1.9. It provides 
the force for all separations between the Casimir-Polder limit for d 2> i?, and the PFA 
result for R/d 1/2. At a surface-to-surface distance 4i?/3 (R/d = 0.3), the PFA 
overestimates the energy by a factor of 10. Including up to Z = 32 partial wave orders 
and extrapolating based on an exponential convergence in I, the Casimir energy has 
been determined down to R/d = 0.49 [32]. The interaction between a sphere and a 
plate has been obtained recently and deviations from the PFA have been quantified [30]. 
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1.3 

Dependence on material properties 



In the previous Sections we liave considered perfectly conducting bodies. For real 
metals with finite, frequency dependent conductivity, or more general dielectric me- 
dia, Casimir interactions are modified. This is a natural consequence of the fact that 
dipole and higher multipole polarizabilities depend on the material properties of a body. 
Thus, the induced fluctuating currents depend not only on shape but also on material 
composition. For practical applications and experimental tests it is therefore important 
to understand the collective effects of shape and material on Casimir interactions. A 
macroscopic theory that fully accounts for the material dependence of the interaction 
between two planar surfaces has been established by Lifshitz [47] in 1956. Until re- 
cently, only approximations, Umited to short separations between bodies or sufficiently 
diluted media, have been available to study the interaction of dielectric media of ar- 
bitrary shapes. The scattering approach described in the previous Section has paved 
the way for studying in detail the material and shape dependence of Casimir forces 
beyond the case of planar surfaces and without the connmon approximations. We first 
provide a simple derivation of the Lifshitz result for two surfaces within the scattering 
approach. Then we focus on an example that is of particular interest to the behavior of 
nano-particles as they appear, i.e., in suspensions where correlations between material 
and shape effects are important. 

1.3.1 

Lifshitz formuia 



Consider two material halfspaces that are bounded by planar, parallel surfaces with 
a vacuum gap of width d between them. The material in the two halfspaces can be 
different and is characterized by the dielectric functions ea(w) and magnetic perme- 
abihties udoj) where a = 1,2 numbers the halfspaces. A compact derivation of 
the Casimir-Lifshitz interaction between the two surfaces follows from the scattering 
formula of Eq. (1.44). The T-matrix of a planar dielectric surface is given by the Fres- 
nel coefficients which are usually expressed in a planar wave basis. When we define 
the two polarizations relative to the surface normal vector, the T-matrix is diagonal in 
polarization and in the wave vector k\\ parallel to the surface. The diagonal matrix 
elements are 

" flc(tCK)p + Pa ^^^^^ 
y _ ta{icn)p - Pa 

(iCK)p -I- Pa 

where p = ^ + fejj and pa = ea{icK)iJba{icK)K^ + fey. The translation matrices 
for translations perpendicular to the surfaces by a distance d are also diagonal in fe|| in 
the planar wave basis and the diagonal elements have the simple form 



,Mfc|| 



Uaa',Ek« = e (1.58) 
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for a a' = 1, 2, i.e., they do not couple E and M polarizations and are identical 
for the two polarizations. The determinant of Eq. (1.44) leads to a product over all fey 
which becomes an integral after taking the logarithm. The resulting Casimir-Lifshitz 
energy has two separate contributions form M and E polarizations (TE and TM modes, 
respectively), 



hcA 



poo pc 

/ dK 
Jo Jo 



kndkn In 



ei {iCK)p - Pl 62 {iCK)p - P2 
ei {iCK)p + Pl 62 {icK)p + P2 



fll {iCK)p - Pl H2{iCK)p-p2 ^-2pd 
111 {iCK)p + Pl IJ.2 {iCK,)p + P2 



(1.59) 



where A is the surface area. This results generalizes the Casimir interaction between 
two perfect metal plates of Eq. (1.5) to dielectric materials. 



1.3.2 

Nanoparticles: quantum size effects 

The Lifshitz formula of Eq. (1 .59), while derived for infinitely extended planar surfaces, 
is commonly applied within a proximity approximation also to curved surfaces of 
particles of finite size. This leads to predictions for the interaction that are limited 
to particles that a very large compared to their separations. To be able to study the 
interaction of particles of arbitrary sizes and separations, a theory is needed that is 
a generalization of the Lifshitz formula to bodies of arbitrary shape. Such a general 
theory provides the scattering formula of Eq. (1.44). The challenges in applying this 
formula consist in the computation of the T-matrix for bodies with general dielectric 
functions and in the proper modelling of the dielectric response of the bodies. The 
latter is especially important for nano-particles for which bulk optical properties are 
modified by finite-size effects. 

Some characteristic effects of the Casimir interaction between nano-particles will be 
discussed in this Section by studying two spheres with finite conductivity in the limit 
where their radius R is much smaller than their separation d. We assume further that R 
is large compared to the inverse Fermi wave vector -ir/kp of the metal. Since typically 
n/kp is of the order of a few Angstrom, this assumption is reasonable even for nano- 
particles. To employ Eq. ( 1 .44), we need the T-matrix of a sphere with general dielectric 
function e(a;) which generalizes the matrix of Eqs. (1.54), (1.55). All elements of this 
matrix are known explicitly, see, e.g., [32]. Relevant to the interaction for d^ R are 
the dipole matrix elements (I = I' = 1) at low frequencies k. To proceed, we need 
information about the dielectric fimction on the imaginary frequency axis lj = icK for 
small K. Theories for the optical properties of small metallic particles [62] suggest a 
Drude hke response 

e(icn) ^ 1 + 4n^^^ , (1.60) 

CK 

where a{icK) is the conductivity which approaches for k ^ the dc conductivity auc- 
For bulk metals adc = uj^t/At: where ojp = •^4e^ft|,/37rme is the plasma frequency 
with electron charge e and electron mass me, and r is the relaxation time. With 
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decreasing dimension of the particle, adc{R) is reduced compared to its bulk value due 
to finite size effects and hence becomes a function of R [62]. 

In the low frequency Umit, with e(icK) of Eq. (1.60) the T-matrix elements for 
magnetic and electric dipole scattering {I = I' = 1) are diagonal in m and have the 
series expansion 

47r^ffdcC5). n fin 

J- MM Imlm = — TF [^KH) +... (.l-Ol) 

45 C 

To leading order ~ the electric dipole matrix elements are identical to those of a 
perfectly conducting sphere and finite conductivity modifies higher orders only. In the 
magnetic dipole matrix elements, however, the leading term — («:i?)^/3 of the perfect 
conductor result of Eq. (1.54) is absent. This is consistent with the observation that 
the magnetic dipole polarizability is reduced by a factor ~ (Kf?)^[e(icK) — 1] and 
e(icK) — 1 ~ due to Eq. (1.60). 

When the matrix elements of Eqs. ( 1 .6 1), ( 1 .62) together with the translation matrices 
Uco,' in spherical coordinates are substituted into Eq. (1.45) an expansion for large 
distance d yields the Casimir energy of two spheres 

E ^ 23fl« ( R<Tdc{R) 45 c \R' 

he 4ir (f \ c 47r2 Rcri^{R) ) d» " " " " ^ ' 

The leading term is material independent but different from that of the perfect metal 
sphere interaction of Eq. (1.56) since only the electric polarization contributes to it. 
At next order, the first and second terms in the parentheses come from magnetic and 
electric dipole fluctuations, respectively. Notice that the term ~ 1 /<f is absent in the 
interaction between perfectly conducting spheres, see Eq. (1.56). The limit of perfect 
conductivity, (Jdc — > oo cannot be taken in Eq. (1.63) since this limit does not commute 
with the low n or large d expansion. 

In order to estimate the effect of finite conductivity and its dependence on the size of 
the nano-particle, we have to employ a theory that can describe the evolution of adc {R) 
with the particle size. A theory for the dielectric function of a cubical metallic particle 
of dimensions R^ ir/kp has been developed within the random phase approximation 
in the limit of low frequencies <g; c/R [62]. In this theory it is further assumed that the 
discreteness of the electronic energy levels, and not the inhomogeneity of the charge 
distribution, is important. This impHes that the particle responds only at the wave vector 
of the incident field which is a rather common approximation for small particles. From 
an electron number-conserving relaxation time approximation the complex dielectric 
function is obtained which yields the size-dependent dc conductivity for a cubic particle 
of volume [62]. It has been shown that the detailed shape of the particle does not 
matter much, and we can set a = {Att/SY^^R which defines the volume equivalent 
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Fig. 1.10 Dimensionless dc conductivity acic(K) in units of 
e^/2?iao (with Bohr radius ao) for a Aluminum sphere with 

ep = 11.63eV, ir/kp = 1.8A and r = 0.8 ■ lO^^-^sec as 
function of the radius R, measured in units of ir/kp, see 
Eq. (1.64). Also shown is the corresponding ratio Ro-^c (fi)/c that 
determines the Casimir interaction of Eq. (1 .63). The bulk dc 
conductivity o-dc(oo) = 17.66 is indicated by the dashed line. 



sphere radius R. This yields the estimate 

kpa/ -K 



inkpa -\- -K^ 487r 



kpa/TT r 

Re m'^{{kFa/n)^ — m^) x < 



with 



7rm / ir 



(1.64) 



—Zmtanzm meven 

+ 2m cot Zm TH odd 



(1.65) 



and where adc(oo) — lj^ r/47r is the bulk Drudedc conductivity and F = {h/TeF){kFa/-K 
is a linewidth with Fermi energy ef- The factor in square parentheses multiplying 
o"dc(oo) describes quantum size effects and leads to a substantial reduction of the dc 
conductivity for nano-scale particles. While the above expression is applicable for 
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Tr/kp -C a, it suggest that for ir/kp — a the particle ceases to be conducting which 
is consistent with a metal-insulator transition due to the localisation of electrons for 
particles with a size of the order of the mean free path. It is instructive to consider the 
size dependence of (Tdc{R) and of the Casimir interaction for a particular choice of ma- 
terial. Following [62], we focus on small Aluminum spheres with ep = 11.63eV and 

r = 0.8- 10~ sec. These parameters correspond to tt/Zcf = 1.8A and a plasma wave- 
length Ap = 79nm. It is useful to introduce the dimensionless conductivity adc{R), 
which is measured in units of with Bohr radius ao, so that the important quan- 

tity of Eq. (1.63) can be written as Radc{R)/c = {a/2){R/ao)adc{R) where a is the 
fine-structure constant. The results following from Eq. (1.64) are shown in Fig. 1.10. 
For example, for a sphere of radius R — lOnm, the dc conductivity is reduced by a fac- 
tor ftj 0. 15 compared to the bulk Drude value. If the radius of the sphere is equal to the 
plasma wavelength Ap, the reduction factor « 0.8. These results show that shape and 
material properties are important for the Casimir interaction between nano-particles. 
Potential applications include the interaction between dilute suspensions of metalUc 
nano-particles. 

1.4 

Casimir force driven nanosystems 

We have seen that Casimir forces increase strongly with decreasing distance and hence it 
can be expected that they are important in devices that are composed of moving elements 
at short separations. Indeed, a common phenomena seen in nano-mechanical devices 
is stiction due to attractive van der Waals and Casimir forces. This effect imposes 
a minimum separation between objects in order to prevent that they stick together. 
However, one can make also good use of Casimir interactions in nano-devices by 
employing them to actuate components of small devices without contact [1, 29, 2]. In 
[29] it has been demonstrated that this can be achieved by coupling two periodically 
structured parallel surfaces by the zero-point fluctuations of the electromagnetic field 
between them. We will consider this effect here as an example for Casimir force 
induced non-Hnear dynamics, providing a direct application of the results obtained 
in Section 1.2.2. We have seen that the broken translation symmetry parallel to the 
surfaces results in a sideways force which has been predicted theoretically [34, 35] and 
observed experimentally between static surfaces [21]. If at least one of the surfaces 
is structured asymmetrically there is an additional breaking of reflection symmetry 
and the surfaces can in principle be set into relative lateral motion in the direction of 
broken symmetry. The energy for this transport has to be pumped into the system by 
external driving. This can be realized by setting the surfaces into relative oscillatory 
motion so that their normal distance is an unbiased periodic function of time. Since the 
sideways Casimir force decays exponentially with the normal distance [see Eq. (1 .41)], 
the surfaces experience an asymmetric periodic potential that varies strongly in time. 

This scenario resembles so-called ratchet systems [58] that have been studied exten- 
sively during the last decade in the context of Brownian particles [41], molecular motors 
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[3] and vortex physics in superconductors [24], to name a few recent examples. Most 
of the works on ratchets consider an external time-dependent driving force acting on 
overdamped degrees of freedom to rectify thermal noise. For nano-systems, however, 
it has been pointed out that inertia terms due to finite mass should not be neglected 
and, actually, can help the ratchets to perform more efficiently than their overdamped 
companions [48]. Finite inertia typically induce in Langevin dynamics deterministic 
chaos that has been shown to be able to mimic the role of noise and hence to generate 
directed transport in the absence of external noise [49]. Here we use this effect in 
the different context of so-called pulsating (or etfectively on-off) ratchets where the 
strengths of the periodic potential varies in time [58]. We consider weak thermal noise 
only to test for stability of the inertia induced transport — not as the source of driving^*. 

It has been demonstrated that the system described above indeed allows for directed 
relative motion of the surfaces due to chaotic dynamics caused by the lateral Casimir 
force [29], The transport velocity is stable across sizeable intervals of the amplitude 
and frequency of surface distance oscillations and damping. The velocity scales linear 
with frequency across these intervals and is almost constant below a critical mean 
distance beyond which it drops sharply. The system exhibits multiple current reversals 
as function of the oscillation amplitude, mean distance and damping. This "Casimir 
ratchet" allows contact-less transmission of motion which is important since traditional 
lubrication is not applicable in nano-devices. This actuation mechanism should be 
compared to other actuation schemes as magnetomotive or capacitive (electrostatic) 
force transmission. The Casimir effect induced actuation has the advantage of working 
also for insulators and does not require any electrical contacts and/or external fields. 
Other applications of zero-point fluctuation induced (van der Waals) interactions to 
nano-devices have been experimentally realized already to construct ultra-low friction 
bearings from multiwall carbon nanotubes [23]. 

In the following, we consider two (on average) parallel metallic surfaces with peri- 
odic, uni-axial corrugations (along the j/i-axis) that have distance H, see inset (a) of 
Fig. 1.11. To begin with, we assume that both surfaces are at rest with a relative lateral 
displacement b. Then the surface profiles can be parametrized as 

oo 

hi{yi) = a^c„e'"'"''i/^i+c.c., (1.66) 

n=l 

oo 

fe(yi) = a^d„e'"™(^^-'')/^^+c.c. , (1.67) 

n=l 

where a is the corrugation amplitude, Ai, A2 are the corrugation wave lengths, c„, d„ 
are Fourier coefficients. 

The dependence of the Casimir energy £ on H and b causes macroscopic forces on 
the surfaces. For a varying separation H this is the normal Casimir attraction between 
metallic surfaces, modified by the corrugations. Below we will assume H = H(t) to 
be a time-dependent distance that is kept at a fixed oscillation by an additional external 

^^In the absence of inertia, finite thermal noise is necessary for on-off ratchets to generate directed 
motion. 
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force from clamping to an oscillator. In such setup, the surfaces can react freely only 
to the lateral force component Tiat{h, H) = —d£/dh. The results of Section 1.2.2 are 
readily extended to periodic profiles of arbitrary shapes as described by Eq. (1.66). The 
corrugation lengths have to be commensurate, Ai /A2 = p/q with integers p, q in order 
to produce a finite lateral force per surface area. For the purpose of this example, it is 
sufficient to consider the case p = 1. Generalizing the result of Eq. (1.37), the lateral 
(6-dependent) part of the Casimir energy per surface area can then be written as 

£{h) = ^f2 (c„d_„,e-^-"''/^^ + c.c.) J (n^\ (1.68) 

n=l ^ ^ 

to order a^. The exact form of the function J(x-) = JTMix) + Jte{x) is given 
by Eqs. (1.38), (1.39). For the present purpose it is sufficient to use the simplified 
expression 

J{x) ~ ^ (1 + 27rx + 7a;' + 32x'') e"^"" (1.69) 

with 7 = 12.4133, which is exact for both asymptotically large and small x and approx- 
imates the exact results with sufficient accuracy for all x. (The maximal deviation from 
the exact result is ^ ±0.5% around x = 0.5.) The Casimir potential of Eq. (1.68) has 
two interesting properties which are useful to the construction of a ratchet. First, it de- 
cays exponentially with H, and thus can be essentially switched on and off periodically 
in time by oscillating H. Second, the potential is not only periodic in b but acquires 
asymmetry from the surface profiles at small H <^ \ and an universal symmetric shape 
forH^X since the effect of higher harmonics of the surface profile is exponentially 
diminished as discussed in Section 1.2.2. 

The relative surface displacement b{t) can be considered as a classical degree of 
freedom with inertia. Its equation of motion is described by Langevin dynamics of the 
form 

pb + jpi) = J^iat[6, H{t)] + V^^m , (1.70) 
where p is the mass per surface area, 7 the friction coefficient, T the intensity (di- 
vided by surface area) of the Gaussian noise ^(t) with zero mean and correlations 
{^{t)^{t')) = 6{t — t') so that the Einstein relation is obeyed. This stochastic term 
describes ambient noise due to effects of temperature and pressure. (Additional contri- 
butions from thermally excited photons to the Casimir force can be neglected at surface 
distances well below the thermal wavelength he/ (2T).) The system is driven by rigid 
oscillations of one surface so that the distance H{t) — Hag{t) oscillates about the mean 
distance Ho with g(t) = 1 — 7] cos(nt). For simplicity, we consider now equal corruga- 
tion lengths Ai = A2 s A. We define the following dimensionless variables: b = b/X, 
s = t/r for lateral lengths and time with the typical time scale r = {X / a) \/~pH^Jhc 
resulting from a balance between inertia and Casimir force. Hence velocities will be 
measured in units of vo = X/t. There are five dimensionless parameters which can 
be varied independently for fixed surface profiles: the damping 7 — tj, the angular 
frequency ui = t51, the driving amplitude rj, the scaled mean distance Ho/X and the 
noise intensity T = (T/Rc)(J?o/o^)- The dimensionless equation of motion for b{s) 
reads 

l + 4 = :Fut[b,g{s)] + ^/2^iis) (1.71) 
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Fig. 1.11 The lateral Casimir force acting between the two 
surfaces as function of the shift b at time s = and half period 
s = tt/llj (drawn to a larger scale by a factor 10^) for parameters 
r] = 0.65, Ho = O.IA. Insets; (a) Surface profiles at their 
equilibrium position at 6 = 0.182 (b) Periodic variation of the 
maximum force at S = with time. 



with the Casimir force 

^i^t{b,g) = ^'^f„cos{2TTnb)J (ng^j , (1.72) 

^ n=l ^ ^ 

where we have chosen surface profiles with c„ = i \//„/(2n), d„ = ^/„/(2n) with 
real coefficients /„ in Eq. (1.66), and ^(s) = 1 — 77 cos((x;s). 

Directed transport is possible in certain parameter ranges even in the deterministic 
case where noise is absent. However, to probe the robustness of transpoit, we consider 
in the following primarily the limit of weak noise by choosing T = 10~^. In fact, it has 
been shown for underdamped ratchets with time-independent potentials and periodic 
driving that even an infinitesimal amount of noise can change the rectification from 
chaotic to stable [48]. To look for similar generic behavior of the pulsating ratchet, 
we consider a specific geometry consisting of a symmetric and a sawtooth-like surface 
profile corresponding to three harmonics with /i = 0.0492, /2 = 0.0241, /a = 0.0059 
and /„ = for 71 > 3. Inset (a) of Fig. 1.11 shows these profiles in their stable position 
with 6 — 0. 182 that minimizes the Casimir energy. The resulting spatial variation of the 
Casimir force with b is plotted in Fig. 1 . 1 1 for minimal (s — 0) and maximal (s — tv/lu) 
surface distance with parameters Ho/\ = 0. 1, 77 = 0.65. It can be clearly seen that the 
asymmetry is reduced at larger distance where the variation of the force becomes more 
sinusoidal. Inset (b) shows the on-off-like time-dependence of the force amplitude at 
6 = due to the oscillating surface distance. 
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Fig. 1.12 Mean {(v)) and standard deviation a of the (negative) 
veiocity as function of the driving amplitude ri for the frequencies 
Ld = 5.0 and = 4.72 (for the latter only the stable plateau is 
shown). The parameters are Ho = O.IA, 7 = 0.9, f = 10~^. 
Inset: Dependence of the same quantities on frequency for fixed 
T] = 0.65. Straight dashed lines correspond in both graphs to the 
velocity uj/{2tt). 



The non-linear equation of motion of Eq. (1.71) has to be solved numerically. The 
trajectory 6(s) was obtained from a second order Runge-Kutta algorithm. As initial 
conditions an equidistant distribution over the interval [—1, 1] for 6(0) and 6(0) = is 
used. For each set of parameters 200 different trajectories are calculated from varying 
initial conditions and noise, each evolving over 4 x 10^ periods 2tv/lli so that transients 
have decayed. The average velocity {{v)) involves two different averages of 6(s): The 
first average is over initial conditions and noise for every time step, then the averaged 
trajectory is averaged over all discrete times of the numerical solution. For an efficient 
directed transport it is not sufficient to have only a finite average {{v)). To exclude 
trajectories with a high number of velocity reversals, the fluctuations about the average 
velocity must be small, i.e., the variance = {{v^)) — {{v))^ must be smaller than 

{{vr- 

Naively, one can expect directed motion of the surface profile h2 {yi ) into the positive 
1/1 -direction {b < 0) since the Casimir force in Fig. 1.11 is asymmetric with negative 
values lasting for longer time than positive ones. However, the actual behavior is more 
complicated due to chaotic dynamics. Fig. 1.12 shows the dependence of the average 
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velocity and its standard deviation a on the driving amplitude rj and frequency ui for 
Ho = O.IA, 7 — 0.9. For a fixed frequency there is an optimal interval of driving 
amplitudes across which the average velocity is almost constant with {{v)) ~ — a>/(27r). 
Small deviations from the latter value result from noise as has been checked by studying 
the dynamics at T = 0. At higher driving amplitudes a second narrower interval with 
maximal {{v}) is observed which is more strongly reduced and smeared out from its 
deterministic value — 2 x w/(27r) by noise. At the plateaus of constant velocity the 
standard deviation a is substantially reduced, rendering transport efficient. Outside the 
plateaus velocity reversals occur and a increases linearly with 77. For fixed amplitude 
77, the average velocity is stable at the value — a;/(27r) over a sizeable frequency range 
(see inset of Fig. 1.12). 

In order to understand the observed behavior it is instructive to analyze the dynamics 
in the three dimensional extended phase space. Attractors of the long-time dynamics 
can be identified from Poincare sections using the period 2it /lo of the surface oscillation 
as stroboscopic time. To obtain a compact section, the trajectory is folded periodically 
in j/i on one period of the Casimir potential. From these sections we can distinguish 
between periodic and chaotic orbits. As a start, we consider the deterministic limit 
with T = 0. The plateaus around rj = 0.65 and rj = 0.7 result both from periodic 
orbits of period one, corresponding to a single point in the Poincare section. On the 
right (downward) edges of the first plateaus we observe period doubling, i.e., a periodic 
attractor with period two. Upon a further increase of 77, chaotic orbits dominate the 
motion. Hence the system exhibits a period-doubling route to chaos with enhanced 
velocity fluctuations. The findings apply basically also to weak noise (f = 10"^) 
but the sharp points of the periodic attractors in the Poincare sections are smeared out 
leading to a decreased ({v)) . The transition from chaotic to periodic dynamics at the 
beginning of the rising edge of the plateaus is accompanied by a velocity reversal. This 
is consistent with the observation for non-pulsating potentials that velocity reversals 
are due to a bifurcation from chaotic to periodic dynamics [49]. 

The amplitude of the Casimir potential can be tuned by varying the mean distance Ho . 
From Fig. 1.13(a) we see that the dynamics show a sharp transition at a critical Ho/\ 
from efficient transport with large {{v)) and small cr to chaotic dynamics with vanishing 
velocity. The transition is accompanied by a velocity reversal and peaked velocity 
fluctuations. Interestingly, below the transition {{v)) is almost constant independently 
of Ho/X. The observed transport behavior is also stable against a change of effective 
damping 7 as shown in Fig. 1.13(b). Whereas fluctuations increase with decreasing 
7, there is a stable plateau of constant average velocity across which fluctuations are 
diminished. In the deterministic limit, additional plateaus with inverted and doubled 
average velocity are observed by varying 7 and rj. Remnants of a second plateau around 
7 = 1.9, washed out by noise, can be seen in Fig. 1.13(b). 

It is interesting to estimate typical velocities vo = A/r. With the typical lengths 
Ho — O.l^m, a = lOnm realized in recent Casimir force measurements [21] and an 
area mass density of p = lOg/m^ for silicon plates with a thickness of a few microns, 
one obtains vo = y/hca?/pH^ w 5.5mm/s. The actual average velocity voijj/2'n is 
of the same order for the frequencies studied above. For A = l/xm, the time scale is 
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Fig. 1.13 Mean ((v)) and standard deviation a of the (negative) 
velocity as function of (a) Vne mean piate distance Ho for 7 = 0.9 
and (b) damping 7 for Hq = 0.1 A. The other parameters are 
V = 0.65, Lo = 5.0, f = 10-3. 

T — \/vo ^ lO^^'s leading to driving frequencies and damping rates in the kHz range 
for the parameters considered here. 

The results show that Casimir interactions can offer novel contact-less translational 
actuation schemes for nanomechanical systems. Similar ratchet-like effects are ex- 
pected between objects of different shapes as, e.g., periodically structured cylinders, 
inducing rotational motion. The use of fluctuation forces appear also promising to 
move nano-sized objects immersed in a liquid where electrostatic actuation is not pos- 
sible. Another application is the separation and detection of particles of differing mass 
adsorbed to the surfaces. For surfaces oscillating at very high frequencies additional 
interesting phenomena related to the dynamical Casimir effect occur [37], leading to 
the emission of photons that could contribute to ratchet-like effects as well. 

1.5 

Conclusion 

There has been recent interest in applying Casimir interactions to the design of nanome- 
chanical devices [17, 16, 10, 9, 11]. In such devices as sensors and actuators, attractive 
Casimir forces can strongly influence their function due to unwanted stiction between 
small elements at nano-scale separations. But one can make also good use of Casimir 
interactions in actuators where they can lead to interesting non-linear dynamics. Re- 
cently also repulsive Casimir forces between bodies in a liquid, predicted some decades 
ago by Lifshitz, Dzyaloshinskii and Pitaevskii for planar surfaces [26], have been mea- 
sured between a sphere and a plane [52], suggesting a way to suppress stiction. Hence 
it is important to understand the dependence of Casimir forces on shape and material 
properties beyond common approximations that apply only to weakly curved surfaces. 
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This conclusion is corroborated by the relevance of Casimir interactions to a plethora 
of phenomena such as wetting, adhesion, friction, and quantum scattering of atoms 
from surfaces. In this Chapter, some characteristic effects of shape and material on 
Casimir interactions have been presented using the examples of geometries that are 
typical to nanosystems. Most of the presented results could be obtained only recently 
by newly developed theoretical tools that have been described here. The important 
study of correlations between shape and material effects and the additional implica- 
tions of interacting fields in Casimir effects due to critical fluctuations [39] are largely 
unexplored. It is expected that the recent progress on the experimental and theoretical 
side will unveil novel phenomena and provide a better understanding of fluctuation 
induced interactions with interesting impUcations for nanosystems. 
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